Theory of real space imaging of Fermi surfaces 
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A scanning tunneling microscope can be used to visualize in real space Fermi surfaces with buried 
impurities far below substrates acting as local probes. A theory describing this feature is developed 
based on the stationary phase approximation. It is demonstrated how a Fermi surface of a material 
acts as a mirror focusing electrons that scatter at hidden impurities. 

Introduction. The scanning tunneling microscope (STM) and spectroscopy provide a unique way of vizualising 
quantum effects from the most basic to the extremely complex ones. Among a wider range of capabilities, STM allows, 
for example, the investigation of electronic wave interferences after electrons scatter at defects on surfaces. Ring-like 
charge oscillations have been observed around Cs adatoms on Ag(lll)Q, electrons waves might even be confined in 
corrals0-0| and used for quantum holography^. These oscillations are important to understand since they mediate, 
for example, interactions between atoms sitting on surface 0, 0|- 

In contrast to surface impurities, research on subsurface defects has been less intense because of the inherent 
experimental and theoretical difficulties involved in the investigation. Recently, a strong stream is being created 
towards the ultimate goal of understanding buried defects and their accompanying electronic waves interferences 
0- 

In particular, we have shown that these interferences can be surprisingly localized and anisotropic on a real 



o 



svstem [14Hl8| |. Using STM and first-principles calculations, we have demonstrated that such effects are induced 
by the shape of the Fermi surface of the bulk substrate, if we probe the Fermi energy in the experiment. Considering 
O i copper as a host and cobalt as impurities, we showed that the very simple Fermi surface of copper bears very flat 
regions that cause, surprisingly, strong anisotropy of the screening charge distribution. For instance, our observations 
allowed us to conclude that Fermi surfaces can be vizualised in real space with STM. Additionally, we proposed to 
utilize a buried impurity surfaces as a local probe for a nanosonar device that is able to map buried defects and 
interfaces and many of their properties, e.g. electronic or magnetic 

This is not unique to copper since nature is full of more complex Fermi surfaces that could lead to unforeseen 
. consequences. Besides the works cited earlier, there have been other works discussing the physics of buried effects. 
Brovko et al. [3] discuss the interesting possibility of detecting the magnetism of large Co-nanostructures buried below 
Cu(lll) surface. With the nanostructure size chosen, though, it is difficult to untie and see a focusing effect and 
observe the bulk Fermi surface. By developing models to calculate the conductance measurable by STM, Avotina 



and co-workers [20| performed a thorough investigation of the effect of Fermi surface shape. Garcia- Vidal et al. [21 1 
and Reuter et a/.[22| have observed similar focusing effects in Au/Si and CoSi2/Si interfaces from another perspective 
• *"H , provided with ballistic electron emission microscopy. To understand their results, they have developed a theory based 
/\ on the Keldysh formalism. 

Our goal is to present a demonstration and a simple theory that takes into account the bandstructure of the host 
and the coupling of the electronic states to the defects they scatter at. Our theory is based on an analysis of Friedel 
oscillations around impurities by using the so-called stationary phase approximation valid for large distances away from 
the impurity. Discussing the effect of the Fermi surface on the electronic propagation goes back to the seminal work 
of Roth and co-workers [23| who developed a theory for the Ruderman-Kittel-Kasuya-Yosida (RKKY) interactions [24j 
mediated by anisotro pic Fermi surfaces. This theory was behind the oscillation behavior of the interlayer magnetic 



exchange coupling[25-27] in terms of calipers of the interlayer Fermi surface. Our aim is to describe the local density 



of states (LDOS) at few Angstroms in the vacuum above the surface that bears a buried single impurity. According 



to the model of Tersoff-Hamann [28[, the LDOS is proportional to the experimental scanning tunneling spectra signal. 
Before discussing the effect of the surface we will develop the theory for a pure bulk material. 

Theory of Asymptotic Behavior. As mentioned previously, our goal is to get the asymptotic behavior of the 
charge oscillations far from the impurtiy. We use Green functions since they allow to work with the Dyson equation. 
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Once known, the space and energy resolved charge can be extracted from the Green function G(r,r';E) for r = f l : 

n(f;E) = --lmG(f,r;E). (1) 
The Dyson equation giving the Green functions of an ideal host perturbed by an impurity can be written as following 

G(?S',E) = G(r,7;E) + [ / dr{dr 2 G{r,r{-E) x 



x t(ri,r 2 -E)G{r 2 ^-E), (2) 

O 

where G is the Green function of the host and t the t- matrix of the impurity being related to the potential perturbation 

o 

AV(r) = V(r) — V(r) induced by the impurity and being given by (in formal notation) 

t = AV 1 (3) 

1-GAV 

We use cell-centered coordinates by replacing f and f by R + r and R! + r 1 , where R and R' are lattice vectors 
and r and f 1 positions in the unit cell. In the following, we assume that the impurity potential AV(r) and the related 
t-matrix t(f,r';E) are localized in the cell of the impurity. The extension to a more extended perturbation is 
straightforward. 

o 

The unperturbed Green function G of the ideal crystal in eq. [5] can be represented by the spectral representation 

^ i r -*j?„(r + R)m (f + r') 

G(f+R,r' + R';E) = '£-L dk kA ^ , (4) 

with Vb being the volume of the Brillouin zone, v the band index and taking the limit e — > +0. Using additionally 
the translation symmetry of the Bloch functions 

*~ kv (?+R) = e rk ^U- kv {f) with U Sv (f+^ = U Sv (f) . (5) 

We can rewrite eq.[2]for the difference Green function AG(r + R, r + R; E) = G(r + R,r + R;E)-G(r + R,r + R;E) 
or AG R (f, f; E), assuming the impurity at position R = 0, as 



AG (f,f;E) = j j df 1 df 2 G(f+R,r 1 - 1 E)t(f u f 2 - 1 E)G(f2,f+R;E) (6) 
Using the spectral representation from eq. [4] we can formulate as 

AG n (r,r;E) = ^ -1 / f dkdk' 



J,k(r+R) e -ik' (r+R) 



(E + ie-E %v )(E + ie-E % , ul ) 
x U^^iE)U* ev Xf) (7) 



with the t-matrix elements f~g (E) given by 

where we integrate over the volume Vq of the unit cell at the impurity site 0. Here (E) is the t-matrix of the 
impurity describing the scattering process at the impurity of an incoming Bloch wave (fc', v') into an outgoing one 
with (fc, v). 
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In order to evaluate the difference Green function AG R of eq. [7] at large distances R away from the impurity at 

° -* ° -* -* 

position R ~ 0, we analyze G(r + R, r\]E) as well as G(r2, r + R) for large distances R 



V B J E + ie-E iv 

To be able to evaluate the integral by the method of stationary phases, we replace the denominator by an integral 
over the time t: 

■ dt ei(E+u _ ESi/H (1Q) 



E + ie - Ej: v J ih 



The resulting double integral over k and t 

G(f+R,r i; E) = ^2 — jdkj^ ^e^(^)^(r)^(fx) (11) 

is dominated by a phase <p^{k, t) varying fastly as a function of k and t 

<t> v (k,t) = kR+(E + ie-E % ji (12) 

Thus for large R, and connected with this are large times t, the factor e 1 ^ oscillates strongly so that important 
contributions to the integral arize only from regions close to stationary points of (f>(k,t) being given by 



-Z=0 = R- T — and -^ = = E-E s 13 
dk ft dk dt 



Thus contributions to the integrals are only expected from fc-points with E^ v values close to the energy E and 

dE- -* -* 

group velocities v% v = |" with directions close to R and further times t close to R — v^ v t. For the evaluation, 
we first devide the k- integral into a two-dimensional integral over the constant energy surface E? — const, and a 

dE^ 

one-dimensional integral dk z perpendicular to this surface, with the direction z given by the gradient |" . For the 
integration over k and t we expand the phase <j> v (k,t) to second order in small deviations Akj and Atj around these 
stationary points kj and tj 



4> v {kj + Akj,tj + Atj) = kjR 

nEflfa,.^^ (14) 



2 h dk a dkp 



Ak a jAkpj 



In the stationary phase approximation the integral ((9]) can be evaluated analytically, provided expansion (|14[) is 
used for the phases 4' v {k 1 t) and the wave functions ^!% v and are replaced by the values at the stationary points 

kj. By determining the stationary points of all bands v the index j includes also the band index in the following 
derivations. In particular, the integration over the time Atj gives, when the integration is extended from — oo to +oo: 

+00 ^ e -» ; A t > t3 = _^i 5 r Ak ) (15) 

-oo ftVj 

for the ^-component of Akj in the direction of R coinciding with the direction of the group velocity Vj = -^5-. For 
the integration over the energy plane perpendicular to Vj we introduce new coordinates Ak X j and Ak y j such that the 
mass tensor gA ® g Ak is diagonal. Using moreover the identity 

dxe- ltix " = j7^e- 4 * SIgntf (16) 
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we obtain for the Green function (fTTj) for very large R values: 

G{r + R,ry,E) - 



^ V B R\ 



d 2 E d 2 E 



dk % dk lj 



(17) 



Here we have to sum over all critical points kj compatible with the direction R. Thus the Green function varies 
asymptotically as ^ (as the free electron Green function) and oscillates with a factor e * sR . In addition, there is a 
phase factor (fj 



7r r . ,a 2 E. . ,a 2 s N1 
•4 {sign fe )+slgn( iT )} 



(18) 



(y5j is respectively equal to and \ when k Z j is a maximum, a saddle point and a minimum of the surface 
constant energy. Most important is that the Green function has an amplitude determined by the inverse square roots 

q2 rp TP 

of the curvatures Jp- and %nr~ at the critical point kj . In the above derivation, we have implicitly assumed that the 

yj 

contribution from the previously mentioned second derivatives does not vanish at the critical point kj which in general 
is realized. However, if vanishing derivatives occur, we speak of a higher order critical point, which results in an even 
slower decrease of the Green function than i. Since such critical points are very important for STM observations, 
leading for subsurface impurities to strong intensity "spots" in certain directions, we will discuss these anomalies in 
the upcoming text. 

For the evaluation of the difference Green function AG(r + R, f+ R; E) of eq. [7]we need the asymptotic expansion 

o o 

for both G(f + R,fi;E) as given by eq. [TBI and the analogous one for G(r2,r + R;E). Note that the latter Green 
function describes the propagation from the cell R to the impurity cell, while the former Green function describes the 

back propagation from the impurity to position R. Therefore in the expression for G(r2,r + R;E), in the equation 
analogous to eq. |H]a factor e~ lkR enters, and the role of r and r 2 (replacing j\) has to be exchanged. Analogously, 

the critical points kji of (f> u {k, t) for G(r*2, r + R; E) are, up to a minus sign, identical with the one for G(f + R, r\]E): 
kj' — —kj. Taking all this together, we obtain for the difference Green function 



AG(r + R,r + R;E) = £ (r) 

33' B 



X tr. 



k 3 k 3 \ E )- R 2 

d 2 E d 2 E d 2 E d 2 E 



J(kzj -k zj i)R+i{ipj +ipji) 



dk2 V3 



dk l3' 



dk V3' 



(19) 



For the evaluation of the energy dependent change of the charge density An(f +R; E) we have to take the imaginary 
part of (|19j) . Assuming only one pair of stationary points j and j', this leads to 



An(r + R; E) 



16tt 



■*c (^)I 2 I^ i I 

sin (2k Z jR + 2(fj +S%.) 



R 2 



d 2 E d 2 E 



(20) 



where Sr, is the phase of the i-matrix tt._t.- In order to get the change An(r + R) of the charge density, we have to 
integrate the strongly varying part over the energy of the occupied states 



dE sin (2k zj R + 2tpj + 8% v ) 



-hv-i 



dk z j sin (2k Z jR + 2<pj + Sr ) 



2R 



cos(2k^fR+2(pj + S n .) 



(21) 
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and can replace the energy in all other parts by Ep. Thus An(r + R) is given by 



. _ 16A J|T 2 cos {2k zj R + 2ipj + 6 % ) 

An(r + R) = __^|^ ( r )|2|^ 



3 2 £ d 2 g 



(22) 



At this stage, we can discuss the physical implications of the previous formula. We see from the last equation 
eq. (1221) that the denominator is a crucial factor. If the denominator is very small meaning that the constant energy 
surface has a flat region, big values of the charge density are obtained from this /c-region, leading to a strong focusing 
of intensity in space region as determined by the group velocity Vj. In other words, the curvature of the constant 
energy surface defines the focusing of the charge. 

Few ab-initio results. Let us now go back to the experimentally investigated case: a Co impurity sitting below 
the surface of Cu(lll)[l3|. We use the Full-potential Korringa-Kohn-Rostoker Green function (KKR-GF) method 
which is ideal for investigating impurity problems in real space exactly in the geometry given by the system of interest. 
For instance, we calculated the space-resolved local density of states (LDOS) in the vacuum region integrated over a 
small energy region around Ep (from Ep — 0.136 eV to Ep) at w 6.1 A above Cu(lll) surface and at ~ 3.5 A above 
Cu(001) surface. Below the Cu(lll) surface we have considered a Co impurity at the 6 th and at the 3 rd layer below 
the surface, and below the Cu(001) we consider an impurity at the 8 th layer below the surface. 



(a) 
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FIG. 1: Fermi surface of Copper represented along three directions: (111) is shown in (a), (001) in (b) and (110) in (c). The 
inverse mass tensor corresponding to the denominator of eq. (|19|) is represented by the color in units of the inverse electron 
mass. Small values represented in red lead to high intensities of the charge variation. 

The Cu Fermi surface is rather spherical apart from the band gaps in the L-directions; flat areas with strongly 
reduced curvature are present in the (110) directions enclosed by the two (lll)-necks and two elevations in the (001) 
directions (see Fig. [1]). These flat regions are represented along the three directions (111), (001) and (110) in red in 
Figffl The colour scale on the Fermi surface represent the strength of the inverse mass tensor (denominator of eq. (|19j) ) 
which measures the flatness. This explains the anisotropic charge ripples observed experimentally and calculated from 
first-principles Interestingly, along the (111) direction, the neck of the Fermi surface defines a forbidden region 
where no electrons can scatter explaining the flat region in the charge density changes at the center of Fig. HJa) and 
Fig. EJb). The latter figures arc different since they are produced by an impurity buried at different distances from the 
surface: in (a) it is at 6 layers below the surface while in (b) it sits much closer to the surface at the 3 rd layer under 
the surface. It is interesting to note that when the atom is closer to the surface, the anisotropy of the oscillations 
seems to loose in intensity which is induced by the stronger scattering of the surface state electrons present on the 
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(111) surface of copper. Those surface state electrons are associated with a nearly isotropic two-dimensional circular 
Fermi surface. In Fig. [2jc), parts of the Fermi surface along the (001) direction (see Fig. [ljb)) is probed. This is 
perfomed by assuming a buried impurity in a Cu(001) sample. To improve the visualization of the curvature change 
observed on the Fermi surface of copper, the Fermi surface is colored with blue and green in Fig. 0|c). Here, the 
two regions with opposite curvatures are obvious and are sepated by a region with a low curvature that induces the 
focusing effect. 



( a ) ^^^^^^^ 




(c) 

V 











FIG. 2: (a) and (b): Impurity induced charge density around Ef at a height of ~ 6.1 A above the Cu(lll) surface with an Co 
impurity sitting in the 6 th layer (a) and in the 3 rd layer (b) below the surface; (c): The case of an impurity buried at 8 layers 
below a Cu(001) surface is shown. Red/blue color means enhancement/reduction of the local density of states at Ef- 

Multiple critical points. While the formula (|18p and (fH?|) are valid for an arbitrary number of critical points, we 
have analyzed in the last section only the results for one pair of critical point. If for a given R- value there exist more 
critical points (fci, k2, ...) on the Fermi surface with group velocities v\, tfj, ... parallel to R, than the Green function at 
Ep as well as the charge density An(f; Ep) exhibit several oscillation periods as a function of R. For instance, for two 
points k\ and k% with v\ and V2 parallel to i?, there are oscillation periods determined by the projections k z \ and k Z 2 
on the direction of R(see Fig.[3ji). Due to the double sum over j, j' in (fl~9|) and (|20|) the charge densities An(r; E) and 
An(r) than show three periods being determined by the k z -vahies 2k z i, 2k Z 2 and k z \ + k Z 2- The amplitude of these 
oscillations are determined by the curvatures at these Appoints as well as the wavefunctions ^.(r) and the t-matrix 

elements (E). It is easy to show that for N critical points the number of periods is ^L^+il Sometimes, e.g. for 

symmetry reasons, some of these periods can be the same, e.g. k% and k 2 can be different, but might have the same 
z-component, so that only one period k z \ = k z2 exists. 

The behavior is illustrated in Fig. [3] for two Fermi surface sections. Fig. [3Ja) shows an ellipsoidal Fermi surface, 
which for each given R vector has only one pair of critical point k with group velocity v//R. Also the two vectors kj 
and kji = —kj are shown, the projection of which on the direction R gives the diameter of the Fermi surface (thick 
line) which determines the oscillation period along the direction R. 

Fig- EKb)-(d) show a more complicated Fermi surface, resembling the "dog bone" of the Cu Fermi surface. For the 
direction R shown in Fig. [3][b), three different kj points with Vj || R exists, leading to a total of 6 different periods 
(2ki z , 2/c2z, 2fc3 Z , k\z + fc2 Z , k\ z + &3 2 , k2 Z + k^ z ) for the charge density. On the other hand, if R is perpendicular to the 
main axis as in Fig. EJc), then k z i = k z % and only three periods exist (indicated by the thick lines) while if R points 
along the main axis, there is only one solution. Thus for a given Fermi surface the situation can be quite complex. 

Higher order critical points. If one of the second derivatives or in eg . [T9l or [20l [22l vanishes, than the 

Green function expression and the charge density diverge, meaning that asymptotically these quantities decrease even 
with a smaller exponent than and ^ respectively. We consider here four such cases, dropping j as index of of 

kj and tji 

a)^r = 0, but 0. In this case, we expand the phase factor of eq. [14]for Ak x up to (Ak x ) 3 : 

1 t f) 3 E At ffi E 1 

(f)(k + Ak.t + At) = k R---—Akl---—Akl--- V Ak z At . (23) 
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FIG. 3: Examples of Fermi contours for illustrative purposes, (a) show an ellipsoidal Fermi surface, which for a given R has 
one critical point with a group velocity v parallel to R. In (a) is shown a thick line which length defines the oscillation period. 
The Fermi surface sections (b), (c) and (d) are similar to the "dog's bone" of copper Fermi surface. For an oblique orientation 
of R with respect to the contour's long axis represented in (b), several critical points contribute to the oscillatory behavior. 
In (c), however, R is parallel to the short axis, which reduces the number of critical points. In (d), red disques show possible 
inflection points corresponding to higher order critical points discussed in the text. One possible direction of R probing this 
region is shown. 
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The integration over Ak x , Ak y , and Ak z , as well as the i-integration can then be performed, leading to a Green 
function: 



G(r + R,f i; E F ) 



1 



1 



(24) 



Thus the decay for larger distances is slightly slower than j^. The charge density at Ep, An(r + R;Ep) varies then 
as -V (instead of -s?), while the total charge An(r + R) varies as -\ instead of the familiar 4y of typical Friedel 

PI'S D T IX 



oscillations. For the case of copper, we observe such a situation as shown in Fig. SJa) and (b): Here the color scale 
on the Fermi surface represents the inverse masse tensors strength along the x and y directions that arc defined as 
diagonal elements of the inverse mass tensor of eq. [T5] The directions x and y are chosen this way that is always 

larger than the y-component. It is interesting to note that in Fig. HJa), is always positive for Cu, whereas ^j§- 

shown in Fig. EJb) changes sign and is equal to zero on the green line. The case considered here corresponds to the 
border line between the green and blue areas in fig. |4jc) where one of the second derivatives changes sign. 






FIG. 4: Fermi surface of Copper showing the diagonal elements of the inverse mass tensor in x and y direction in (a) and 
(b), respectively. In panel (c) the Fermi surface of copper is colored following the curvature of the Cu-Fermi surface, which 
determines the amplitude and the phase of the oscillations in eq. [18] and eq. 1201 Blue indicates an inward curvature like the 
curvature of an ellipsoid, green means an outward curvature. At the boundary line between the green and blue areas the 
Gaussian curvature vanishes, indicating higher order critical behaviour as in Eq. 1241 

b) In case that ^0- = = 0, the charge density decreases again slower. Assuming that the third order derivatives 
have radial symmetry, the difference of the Green function varies proportional to — |-, An(R; Ep) decays as — r, and 
the total charge density change drops as -V. Such a case could be induced by the inflection points represented as a 
red disk spot in Fig. EJd). 

c) We consider now a case where the energy E(k x ,k y ) is constant along a k y -\me of length I perpendicular to R, 
with constant group velocity pointing along R. Then AG oc - — ; 1 - i • / and An(R; Ep) oc ~, An(R) oc j^. 







(r\ ° 2e 1 


>* 







The decrease corresponds exactly to Friedel oscillations in two dimensions {x and z), the third dimension k y gives 
just a constant factor I, respectively I 2 for the charge density. Thus the decrease is very slow. 

d) Finally, we consider an extreme case, that E{k xi k y ) is constant on a whole plane with edge length I representing 

o 

a perfectly flat part of the Fermi surface. Then the k x and k y integrations give just a factor of / each, such that G oc I 
and An(R; Ep) oc I and An(R) oc ^. Thus the Green function AG and the charge density have a purely oscillatory 
behavior with no decrease in R, while the Friedel oscillations decrease only as This behavior is typical for a one 
dimensional system, and represents the slowest possible decrease. Most important is that the amplitude varies with 



9 



fourth power of the length I. Thus planar pieces of the Fermi surface lead to particular slowly decreasing oscillations, 
and this effect increases strongly with the size of the platelet. 

Role of the surface. Here we would like to comment on the role of the surface. Up to now we have assumed a 
bulk system. However, as mentioned previously, the quantity of interest for the interpretation of STM experiments 
is the LDOS inside the vacuum at a certain distance z from the surface. According to the theory of Tersoff and 
Harmann[28| this distance corresponds to a position in the center of the STM tip so that tunneling parameters (bias 
Voltage, setpoint current) as well as the tip radius have an impact on z. There are different ways to implement the 

o 

surface into the model. For simplicity, we will work with the spectral representation of the host Green function G 
from eq. (U]). 

The tunneling current has contributions of multiple states, which decay differently into the vacuum. If we assume 
full translational invariance parallel to the surface, the parallel component of wave vector fc|| is conserved. Inside the 
vacuum the electrons obey the Schrodingcr equation of a free-electron: 

|-(fcjj -k 2 ) = AE-<P (25) 

Here $ is the work function of the material and AE = E — fi the electron energy with respect to the chemical potential. 
From this a fcii -dependent decay constant k can be derived: 



< k \\) = \/j^(*- AE )+ k \\ ( 26 ) 

The above expression shows that states with high /c|| -values have a higher k and therefore decay faster into the vacuum. 
Consequently the STM is more sensitive to states near the center of the surface Brillouin-zone with small ku , while 
short-wavelength contributions to the LDOS oscillations, corresponding to large ku , are suppressed. 
In the vacuum, the wave-function of state is then: 

V(f\\,z) oc exp (ik\\f\\j ■ exp (-z^k^ + K%j , (27) 

where we defined kq = K(k» = 0) = •y/2m($ — AE)/h. In order to get an idea how this effect influences the observed 
patterns we perform a Taylor approximation of eq. [27] up to first order in kfr. 



k 
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This approximation is valid for ku <C kq- For the case of copper this is a good treatment of states having ku < 1.1 A 1 
while states having higher fc|| values are over-suppressed . Inserting |2"51 in |2"T1 gives a simple expression: 

( k fA 

exp (— kz) « exp (— kqz) ■ exp . (29) 

\ 2k J 

This is, apart from a general attenuation of the wave functions amplitude by exp (— Kqz), a Gaussian with a width of 
Ofc = kq/z. This implies that one can relate the wave functions from a smaller distance Z\ to greater distances Z2 by 
a convolution (symbol *) with a Gaussian of width a 2 z = kq/(z2 — Z\): 

*(f[|,z 2 )=exp(- Ko (z 2 -z 1 ))-*(r||,z 1 )*exp^i^ . (30) 

This is a very helpful expression as the effect of the tip-sample distance z can be described by a simple Gaussian 
filtering. If we increase z either by choosing tunneling conditions, where the tip is at a larger distance from the surface 
or by using a blunter tip, the wave functions probed by the tip will be increasingly smeared out. If the approximation 
in eq. [29] is not valid, the convolution has to be performed by a Fourier transform of the last term in eq. [27] but 
the whole effect can still be understood as a kind of smoothing filter. This convolution can also be applied to any 
superposition of wave functions as well as to the Green functions which would then decay similarly. 

In other words, one has to be carefull since additionally to the host Fermi surface the vacuum tunneling modifies 
and preselect interferences. Interferences created in the bulk can be different from those measured with the STM 
above the surface of the material. 
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Conclusion. To conclude, as the Fermi-surfaces of most materials deviate strongly from a spherical shape, the 
corresponding propagation of electrons is anisotropic and could reveal new effects in different materials present in 
nature. For instance, the combination of buried impurities and a scanning tunneling microscope could be used as a 
nanosonar to investigate the interior of materials. We developed and presented a theory behind the focusing effect of 
electronic wave oscillations. Additionally, the effect of different kind of Fermi contours' critical points are discussed 
and the consequences for the decay of charge density oscillations are highlighted. 
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